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We have developed a non-perturbative functional renormalization group approach for the random 
field O(N) model (RFO(N)M) that allows us to investigate the ordering transition in any dimension 
and for any value of N including the Ising case. We show that the failure of dimensional reduction 
and standard perturbation theory is due to the non-analytic nature of the zero-temperature fixed 
point controlling the critical behavior, non-analycity which is associated with the existence of many 
metastable states. We find that this non-analycity leads to critical exponents differing from the 
dimensional reduction prediction only below a critical dimension d c (N) < 6, with d c (N = 1) > 3. 

PACS numbers: 75.10.Nr, 11. lO.Hi 
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After decades of intensive investigation, the nature 
and properties of the phase transition in systems with 
quenched disorder remains a much debated topic. Among 
the controversial issues are the critical behavior and 
the phase diagram of the random field Ising model 
(RFIM) 0. For the RFIM, the diagnosis is well estab- 
lished: standard perturbation theory "exactly" (i.e. to 
all orders) predicts an equivalence between the critical 
behavior of the RFIM in dimension d and that of the 
corresponding pure Ising model in dimension d — 2 ( "di- 
mensional reduction" ) [2, El , which contradicts rigorous 
results showing that the lower critical dimension of the 
RFIM is 2 and not 3 . The failure of dimensional re- 
duction (DR) and conventional perturbation theory is at- 
tributed to the existence of a "complex energy landscape" 
characterized by many metastable states [fjj. However, 
limited progress has been made for developing a theory 
able to cope with metastability in random field systems, 
especially for the physically relevant case of the RFIM 
in d = 3. The most significant step has been the for- 
mulation of a perturbative renormalization group (RG) 
approach for the RFO(N>2)M near the lower critical di- 
mension d — 4 6]. Already at leading order, the calcula- 
tion includes an infinite set of marginal operators near the 
zero-temperature fixed point controlling the critical be- 
havior, hence the name functional renormalization group 
(FRG). This perturbative FRG has been further devel- 
oped to study equilibrium and out-of-equilibrium proper- 
ties of elastic manifolds in disordered medias |^I^ ITolfll| . 
and it has been shown that the large scale behavior of 
these systems is controlled by non-analytic renormalized 
actions, with non-analycities encoding the effect of the 
many metastable states at zero temperature. 

In this letter, we develop a non-perturbative FRG for 
the RFO(N)M by combining the ideas of the perturba- 
tive FRG for disordered systems with the formalism of 
the non-perturbative RG for the effective average action, 
based on an exact equation (ERGE) We imple- 

ment a tractable approximation scheme that allows us 
to recover the perturbative results in the relevant limits 



(most importantly, the perturbative FRG at first order 
in e = cl — 4 for N > 2 0, Q) and to describe, approxi- 
mately but non-perturbatively, the ordering transition in 
the whole d — N diagram, including N = 1 and d = 3. 
This provides a comprehensive picture of the critical be- 
havior of random field systems. We find in particular 
that this latter is controlled by a zero-temperature fixed 
point at which the renormalized effective average action 
is non-analytic (albeit in a more complex way than in the 
random elastic manifold case however, the DR 

prediction for the critical exponents breaks down only be- 
low a critical dimension d c (N) < 6, with d c (N = 1) > 3. 

We start with the standard effective hamiltonian for 
the RFO(N)M in d dimensions with an TV-component 
field x( x ) coupled to a random field h(x) with zero mean 
and variance h^(x)h u (ij) = AS fiU 6(x — y). For conve- 
nience, we derive the ERGE for disorder averaged func- 
tions within the replica formalism, but it could similarly 
be obtained by using e.g. the dynamic formulation. We 
thus consider the "replicated" action 



S n [{Xa},{Ja}] = J d "4^E 



dXa ■ dXa 



2 u 4 



(i) 



2^2 y / Xa-Xb y t Xa-Ja 



a,b=l 



where we have introduced sources acting separately on 
each replica, which therefore explicitly breaks the per- 
mutation symmetry between the replicas. One can 
associate to the above action the thermodynamic po- 
tential W n [{J a }} = \ogZ n [{J a }\ where Z n [{J a }} = 
JI\a=i T> Xa,exp(-S n [{Xa},{ J a,}}), and its Legendre 
transform, the effective action r n [{</> a }] = — W„[{J a }] + 
J d d x 2a=i ^ a ' ^ a w hich is the generating functional 
of the vertex functions. In the following, we drop the 
subscript n. 

To investigate the phase diagram and critical behavior 
of the model, we use an ERGE for the effective average 
action 12]. An effective average action Tk at the run- 
ning scale k is obtained by integrating out fluctuations 
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with momenta q> k via the introduction of an infra-red 
(IR) cutoff function R k (q); Tk continuously interpolates 
between the bare action, eq. JIJ, & t the microscopic scale 
k = A and the usual effective action when k — > 0. ft 
follows an exact flow equation, 

d t r k {{ct> a }} = l/2TTd t R k {q){T^ ) [{<t> a },q] + tR k (q))- 1 

(2) 

where d t is a derivative with respect to t = ln(fc/A), 1 
is the unit matrix with elements (2-K) d 5^ d \q — q^S^S^, 

(2) 

Tj, is the tensor formed by the second functional deriva- 
tive of T k with respect to the fields 4»%(q) and <f>%(q'), and 
the trace involves an integration over momenta as well as 
a sum over replica indices and ./V-vector coordinates. 

Eq. J3J) is a complicated functional integro-differential 
equation that cannot be solved exactly but provides a 
convenient starting point for non-perturbative approxi- 
mation schemes. One such scheme that efficiently deals 
with the momentum dependence of the vertex functions 
is the derivative expansion [l^ . However, disordered sys- 
tems require more because inversion of the matrix involv- 
ing in eq J2J) is a difficult task as far as the replica 
indices are concerned for non-integer values of n. As in 
the perturbative FRG approach [10( , we follow the route 
that consists in expanding all functions of the replica 
fields {4>a\ in increasing number of free replica sums. Il- 
lustrated for the potential Uk({<fi a }) at the running scale 
k (i.e. the effective average action for uniform fields), 
this gives: 



n ^ n 

U k {{<t>a}) =Y, U ^ < t>») ~ o E M<t>a,4>b) 



(3) 



a=l 



a,b=l 



where U k , Vk, • • • are continuous functions of their ar- 
guments and satisfy the replica permutation symmetry. 
When all replica fields are equal, each free replica sum 
brings a factor of n and the procedure amounts to an 
expansion in powers of n (with n — > 0) . 

Within this framework, the simplest truncation for Tk 
that already contains the main ingredients for a non- 
perturbative approach of the RFO(N)M is the following: 



Fk{{<f>a}} = J d d xl ^2 (^Z mjk dct> a ■ dcf> a 

\, a—1 

1 n 1 

a, 6=1 J 



(4) 



with one single wave function renormalization for all 
fields, Z m< k, which is defined as the derivative w.r.t. q 2 
of T k evaluated at q = for a (non-zero) field configu- 
ration </> m that minimizes the 1-replica potential Uk- this 
is the so-called pseudo first-order of the derivative ex- 
pansion With the above truncation that keeps only 
the first two terms of the expansion in the free replica 



sums, the ERGE for Tk, eq. (JSJ, can be reduced to cou- 
pled partial differential equations for the functions Uk(4>) 
and 14(0, 0'), whereas a running anomalous dimension 
is defined as r\ k = —dt log Z^ k . The details will be given 
elsewhere. 

To study the critical behavior associated with the or- 
dering transition and search for fixed points (FP) of the 
flow equations, we introduce as usual renormalized di- 
mensionless quantities. However, anticipatin g th at the 
putative FP is expected at zero temperature [1J3, |l4j , it 
is convenient to make explicit the flow of a running tem- 
perature and the associated exponent. For simplicity, let 
discuss first the RFIM. We define a renormalized disor- 
der correlation function Ak(4>,<p') = d^d^i V k (^>, <f)') and 
a renormalized disorder strength A. mik = A.j~((p m , <p m ) ■ 
A running temperature can now be defined by Tk — 
Z mifc fc 2 A/(A 2 A mifc ): when k = A, it reduces to T A = T 
(since from eq. (JTJ, Z A = l/T and A m , A = A/T 2 ). 
An associated running exponent is obtained from 6 k = 
dtlogTk- By using the definition of r/k, one may alter- 
natively introduce an exponent fjk = —6k + 2 + r]k and 
compute it from the equation fjk — 2r\k = <9(A m <k . Di- 
mensionless quantities (noted by lower cases) appropri- 
ate for looking for a zero temperature FP are then: ip = 
{T k Z k k-( d -Vy/ 2 <i>, u k (cp) = T k k- d U k ((f>), v k {<p,<f/) = 
T 2 k- d V k {<P,<P'), and 5 k (ip,ip') = d^v^,^)- The 
procedure can be extended to the RFO(N)M. It is how- 
ever more convenient in this case to introduce the vari- 
ables p = (f 2 /2 and z = if ■ tp'/ (App') 1 / 2 . In scaled form, 
the flow equations for u k {p) and Vk{p, p 1 , z) read (for sim- 
plicity, we drop the subscript k for all quantities but T k ): 

d t u = (2 - d+j] — fj)u + (d — 4 + fj) pu p 

+ 2v d [(N - l)<Mi + S L lf(w L )] (5) 
+ 2v d T k [(N - l)l d (w T ) + l d (w L )] 

d t v = (4 - d + 2i) - 2fj)v + (d~4 + fj)(pv p + p'v p >) 



{(N 



-1) 



(2pv p - zv z ) S T l d (w T )/p + [2p'v p : 

- zv z )5' T l d (w' T )/p' + v 2 z l d hl (w T , w T r )/(2pp')] + 

(1 - z 2 ) [v zz (5 T l d {w T )/p + 5' T l d (w T ')/p') + 2pv% 
If^wr', w L )/p' + 2p'v 2 p ,Jl 1 (w T , w L ')/p - (v 2 z 
+ 2zv z v zz - (1 - z 2 )v 2 zz )li 1 {w T ,WT l )/{2pp l ) 

+ 2 5 L (v p + 2pVp P )l d (w L ) + 6' L (v p > + 2p'v plp ,) 

liM + Wv 2 pAi( w ^ w l')}} -TkVd 
| (AT — 1) [(2pw p - zv z )l d {w T )/p + (2ft V? - zv z ) 
lf(w' T )/p'] + (1 - z 2 )v zz [li(w T )/p + lt(w T ')/p' - 
(v p + 2pv P p) l d (w L ) + ( V + 2p'vp, p ,)l d (w L ') } 



(6) 
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where the indices p, p' , z indicate derivatives w.r.t. 
p,p',z, and v^ 1 = 2 rf +V/ 2 r(d/2); w t ,l = w t ,l{p), 
w' T L = wt,l{p')i and similarly for St.l and S' T L , 
where w T {p) = u p {p) and w L {p) = u p {p) + 2pu pp (p) 
are the transverse and longitudinal masses, whereas 
Mp) = l /( 2 P) v z(p,P,l) and S L (p) = 2pv pp , (p, p, 1) 
are the transverse and longitudinal disorder correlation 
functions when ip = ip' . For brevity we have omit- 
ted the arguments p and p,p',z of all functions. Fi- 
nally, lp(w) and l pp ,(w,w') are the so-called dimension- 
less threshold functions, that essentially encode the non- 
perturbative effects beyond the standard one-loop ap- 
proximation: their definition and properties are discussed 
at length in ref. 0] . For the present work, we choose for 
IR cutoff function R k (q) = Z m>k (k 2 -q 2 )Q{k 2 -q 2 ) where 
is the Heaviside function. When N = 1 and z = ±1, 
one recovers the RFIM (p and p' being used in place of 
ip and ip'). 

The above flow equations are supplemented by equa- 
tions for r]k and fjk- For lack of space, we just give here 
the equation obtained for 2r]k — fjk = — dt &m,k m the case 
of the RFIM: 

2% - fj k = 2v d {lt{u'l)u'f - 4Zg «)<'C 

o /// i si 2 

+ li(u"){8" + -8' 2 - %8 m - -E m )+lf(u")%- 

1 II 1 " 1 
- T k [l d 2 (u'XS' m - li(u>:)C-5>: - ^8' m + ±£ m )]} 

with = lim^^y^^ - 9 v /) 2 ((5(^, </?')- <y?)) 2 and 
S(<y9) = lim v _ +v /(9 v — d v i) 2 8{ip,ip'); the subscript m in- 
dicates that the functions are evaluated for fields equal 
to (p m and primes indicate derivatives w.r.t. ip. Note the 
appearance of the "anomalous" terms E m and 71- E m that 
can only differ from zero when a non-analycity (a "cusp" ) 
in ip — ip' appears in the renormalized disorder function 
8(ip,ip') when ip' — > (p. Actually, if 5(<p,ip') is analytic 
when ip' — > ip, the flow equations for u(ip) and S(ip,ip) 
can be closed by taking from the beginning the replica 
symmetric limit: this analytic behavior in the vicinity of 
replica symmetry is precisely what is implied by the stan- 
dard perturbation theory. In our formalism, breakdown 
of DR thus implies the emergence of a non-analyticity in 
the renormalized disorder correlation function. Note also 
that if a FP is found to eqs. lloKit . and provided that 
^ = 2 — 77 + 77 > 0, it is at zero temperature and tem- 
perature is irrelevant (albeit dangerously) : indeed, in the 
vicinity of the FP, T k flows to zero as k e when k — > 0. In 
most of the following, we will consider directly the T = 
limit, which allows to drop all terms proportional to T k 
in the above equations. 

Because of their structure, the above non-perturbative 
FRG equations reproduce all perturbative one-loop re- 
sults in their range of validity, in particular the e = 
6 — d expansion at first order and the N = 00 limit; 
a stronger property is that one also recovers the per- 



turbative FRG equation at first order in e = d — 4 for 
the RFO(N>2)M 2 6j; In this case, d = 4 being the 
lower critical dimension, the FP occurs (as for the pure 
system) for a value of p m that goes to infinity as 1/e. 
One can thus organize a systematic expansion in pow- 
ers of l/p m . The longitudinal mass becomes very large 
around p m , and by using the known asymptotic proper- 
ties of the threshold functions for large arguments |12j . 
one can derive the flow equations for p m (obtained from 
eq. (JjJJ and the condition u'(p m ) = 0) and for the function 
R{z) = v(p m , p m , z)/{2p 2 a ) (recall that since 8t,„, = 1 by 
construction, R'(l) — l/p m ). This latter reads 

d t R(z) = (e + 2i])R(z) — C4/2 j(7V — I) [R'(z) 2 
+2R'(l)(2R(z) - zR'(z))] + (I - z 2 ) [ - R'(z) 2 (8) 
+ (I - z 2 )R"{z) 2 + 2(R'{1) - zR'(z))R"(z)] } 

where C4 = 2«4Zf(0) = (167T 2 ) -1 , irrespective of the 
choice of the IR cutoff function; in addition, the expo- 
nent r\ is given by 77 ~ C^R![X). The above equations are 
identical to the FRG equations at order e first derived by 
D. Fisher 0. 

The main advantage of the non-perturbative FRG that 
we have developed is that the mechanism by which DR 
and conventional perturbation theory break down can be 
studied in the whole d — N diagram. Although we have 
not yet obtained the full numerical solutions to eqs. (|5I6|I . 
partial solutions and analyses lead us to propose the fol- 
lowing picture. (I) Except when N = 00 and d > 6, the 
analytic FP found in perturbation theory is never stable 
(more precisely once unstable). (2) The stable FP is char- 
acterized by a renormalized disorder correlation function, 
involving 2 fields ip and ip' , which is non- analytic near 
ip = ip' ~ ip m : for instance, in the RFIM S(ip, ip') is non- 
analytic in (ip — ip') and for the RFO(N)M near d = 4, 
R(z) is non analytic in (1 — z); more generally, the non- 
analycity appears in the variable (ip—ip') 2 . (3) The power 
exponent characterizing the non-analycity increases dis- 
continuously with N and d: e.g., near d = 4 for N > 2 
there is a cusp R'(z) - (f-z) 1 / 2 for 2 < N < 18, a "sub- 
cusp" R'{z) - (I - zfl 2 for 18 < N < 18.045, and so 
on, and at large N the exponent varies as N/2; a similar 
trend is observed with increasing d, and in d = 6 — e, the 
exponent goes as 1/e 2 (a result quite similar to that of 
Feldman yy). (4) Only when a cusp is present do the crit- 
ical exponents change from their DR value, so that there 
is a critical line d c (N) (or N c (dj) separating the two re- 
gions of the d— N plane in which the exponents are equal 
(d > d c (N)) or not (d < d c (d)) to the DR predictions. 
To locate this line, it is sufficient to study the flow of the 
second derivative d 2 v/d((ip — if') 2 ) 2 when ip = if' under 
the assumption of analytic behavior: the appearance of 
a cusp, or equivalently of a term in \ip — ip'\ 3 in v, is then 
signalled by a divergence in the second derivative. To 
simplify the analysis, we have used an expansion in pow- 
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FIG. 1: Location of the critical line d c (N) below which a cusp 
appears and exponents differ from their DR values. 

ers of the fields (including all terms up to 4> 4 ) around (p m : 
u = u 2 (p-p m ) 2 , v = 2vi(y/pf7z- p m ) + v 2 (p+ p' -2p m ) 2 + 
vz{p-p l ) 2 +v i {^p~^z-p la ) 2 +v 5 {^p~p^z-p ra ){p+p' -2p m ), 
where V\ = 1 by construction. The result is shown in 
fig. when d — ► 4, one recovers that the critical N is 
18 (see above and ref. 6J) and for the RFIM one finds 
d c (N) ~ 5.1 (within the present approximation). (5) 
Below d c (N), the T — "cuspy" FP is associated with 
critical exponents differing from their DR value; the sim- 
plest illustration is for the RFO(N>2)M at first order in 
e = d— 4: see fig.|5]where we display r\ and fj normalized 
by their DR value fj = T) = e/(N — 2) as a function of 
N (for N = 3, 4, 5, they agree with those of ref. 0). (6) 
For any finite temperature T/., the cusp is rounded: this 
can be inferred, e.g. , from the term Jfe£ m in eq. Q that 
must stay finite as one approaches the FP. As in the ran- 
dom elastic manifold problem 0, ^| , temperature pre- 
vents the flow from being non-analytic at any finite scale 
and one may expect that the rounding of the cusp in- 
volves a boundary layer as Tk — ► and one approaches 
the FP. 

The present non-perturbative FRG approach of ran- 
dom field systems provides a consistent and global pic- 
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FIG. 2: Exponents n and fj (normalized to their DR expres- 
sion) vs TV for the RFO(N>2)M at first order in e = d - 4. 



ture of the critical behavior associated with the ferromag- 
netic ordering transition. The failure of DR and standard 
perturbation theory comes from the existence of many 
metastable states, but the mechanism by which this oc- 
curs is rather subtle: metastability results from an inter- 
play between ferromagnetic ordering and disorder, and it 
plays a role at large scale because the fixed point occurs 
at T = 0. This interplay leads to an effective renor- 
malized random potential (beyond the bare random field 
term) that displays many minima. The renormalized dis- 
order correlation function, which is the second derivative 
of the second cumulant of this random potential, acquires 
a non-analycity (a cusp in low enough d) as it flows to 
the T — fixed point. The physics of such cusps has 
been discussed in the context of random elastic mani- 

folds BSESEi; but in the present case the cusp only 

occurs when the system is in the vicinity of the minimum 
of the non-random potential. Work is in progress to ob- 
tain the full numerical solution of the flow equations and 
compute the critical exponents, as well as investigate the 
connection of the present theory with replica symmetry 
breaking approaches |l6t Ht| and possible formation of 
bound states |l7lll8j . 

LPTL is UMR 7600 at the CNRS. 
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